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Abstract 

Generalized Langevin equations (GLE) with multiplicative white Poisson noise pose the usual 
prescription dilemma leading to different evolution equations (master equations) for the probability 
distribution. Contrary to the case of multiplicative gaussian white noise, the Stratonovich prescrip- 
tion does not correspond to the well known mid-point (or any other intermediate) prescription. 
By introducing an inertial term in the GLE we show that the Ito and Stratonovich prescriptions 
naturally arise depending on two time scales, the one induced by the inertial term and the other 
determined by the jump event. We also show that when the multiplicative noise is linear in the 
random variable one prescription can be made equivalent to the other by a suitable transformation 
in the jump probability distribution. We apply these results to a recently proposed stochastic 
model describing the dynamics of primary soil salinization, in which the salt mass balance within 
the soil root zone requires the analysis of different prescriptions arising from the resulting stochastic 
differential equation forced by multiplicative white Poisson noise whose features are tailored to the 
characters of the daily precipitation. A method is finally suggested to infer the most appropriate 
prescription from the data. 
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I. INTRODUCTION 



Intense and concentrated state-dependent forcing events may often be modeled as multi- 
plicative random jumps, taking place according to an underlying point process. Unlike the 
additive case, which counts a relatively vast literature [l-6|, state-dependent jumps have 
been less investigated 7Hll|, and usually the state-dependency is assumed to be in the fre- 
quency of the jump occurrence, rather than in its amplitude. The generalized Langevin 
equation (GLE) for white multiplicative noise C(^)) which can be either Gaussian or non 
Gaussian, 

x{t) = a{x,t) + &(x)C(t), (1) 
is ill-defined unless a prescription for the evaluation of the stochastic term b(x)C(t) is speci- 



fied 



12| . While this issue is well understood for Gaussian white noise (GWN) [l3|, a precise 



characterization of the noise prescriptions and a clear connection between the different in- 
terpretations are still missing for other kind of noises. 

The last term in Eq. ([1]) for white Poisson (WP) process can be written as ({t) = 
^p{i^, t) = Wi6(t — ti), where the are the times at which jumps occur, 6 is the Dirac delta 
function, and the probability that n jumps occur during a time interval At is given by the 
Poisson distribution -P„(t) = exp (— z/At)(z/At)"/n!. The jumps heights w are independent 
and identically distributed random variables with a probability distribution function (PDF) 
p{w). We note that the multiphcative case b{x)({t) of Eq. ([1]) is a special case, in which 
the X dependence of a more general state dependent white noise Ci^yt) can be factorized 
out. Note that, while it is always possible to reduce the state dependent noise as in ([T]) for 
GWN, because GWN is fully characterized by its mean and variance, this is not the case 
for the WP process. 

The paper is organized as follows. First we show how different prescriptions corresponding 
to the Ito (I) and Stratonovich (S) interpretation of a stochastic differential equation (SDE) 
arise naturally for multiplicative jumps, depending on the relevant time scales of the process. 
In section 3 we present the Master Equation (ME) for a GLE with multiplicative compound 
Poisson process in both the I and S prescriptions. The core of this work is presented in section 
4 where we show how in the linear case, b{x) oc x, the difference between prescriptions 
is properly interpreted as a transformation of the jump size PDFs. We demonstrate the 
relevance of these effects on a minimalist model of soil salinization, describing possible long- 
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term accumulation of salt in soils in arid and semi-arid regions. In this problem, the salt 
mass balance equation is characterized by state dependent losses concentrated in negative 
jumps due to the leaching of salt produced by intense rainfall events. The stochastic equation 
is solved analytically obtaining explicitly the jump distributions that arise in connection to 
the different noise interpretations. 

II. CONNECTION BETWEEN DIFFERENT PRESCRIPTIONS OF A CLE AND 
TIME SCALES OF THE PROCESS. 

We begin with a pedagogical example of a particle that experiences multiplicative im- 
pulsive forcing events, proportional to Or(i), of duration r, in a field characterized by a 
friction coefficient ip. Our analysis is inspired by the work in references [isl. We choose 
e^(t) = d{t/T) with d{z) -> 1 (0) in the limit z oo (-00) so that e^(t) 5{t) in the 
r — i- limit (in the distribution sense). We first consider the case of a single jump event at 
t = to > 0, where the dynamics is described by the Newton equation 

mx{t) = —■ipx+ ipb{x)wQrit — to), (2) 

where the random jump w is drawn from the jump size PDF p{w). Thus in Eq. ([2]) we have 
two time scales, a = m/ip and r. The former is associated with the relaxation time toward 
stationarity, while the latter is related to the characteristic duration of the impulsive forcing. 
Different prescriptions of Eq. ([1]) arise depending on how the two emerging timescales a and 
r in Eq. ([2]) go to zero, i.e. a — t- followed by r — )■ or viceversa (see Fig. [1]). For this 
reason, writing x{t) = b{x)w6{t — to) is ambiguous, being the result of two different limit 
procedures with different physical and mathematical meaning. 

When (7 ^ r and then the zero limit of r is taken in Eq. ([2]), the S p rescription of 



the SDE (11]), which preserves the usual rules of calculus, is obtained [16|, ll7|. For example 
if b{x) = X, the resulting S-equation dln{x)/dt = wQr{t — to), after performing the limit 
r — )• 0, has formal solution x(t) = (1 + 0(t — to){e^ — 1)) xo, where xo = x{0) and 6 is the 
Heaviside function. The corresponding PDF is 

p (H-)) T 

/(x, t) = 5ix - xo)(i - e(t - to)) + e(t - to) ^ '° ^ e(-), (3) 

X Xo 

with initial condition p^{x, 0) = 6{x—xo). If otherwise r ^ cr, then Eq. ([2]) becomes ax+x = 
b{x)w6{t — to). Imposing the conditions of continuity and right and left differentiability in 
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FIG. 1: (Color online) Comparison between trajectories of a particle that undergoes impulsive 
multiplicative forcing in a viscosity field for different timescales (r and a), and the trajectories 
that result from the SDE x{t) = —x{t) Yli=i'^i^i^ ~ interpreted in the I and S prescriptions. 
The jumps in this case are given by Wi = ±0.4. 

to, the initial conditions a;(tQ ) and x(to ), and taking the limit cr — )■ 0, the solution is (again 
for the case b{x) = x) x(t) = xq + XowQ(t — to)- Note that the latter corresponds to the 
solution in the Ito prescription of the SDE (j2]). From the formal Ito solution of Eq. (|2]) we 
obtain the corresponding PDF in the I sense 



p^x, t) = 5{x - Xo)(l - e(t - to)) + 0(i - to 



P 



X — Xq 



Xo 



(4) 



The latter equation can not be made to correspond to Eq. ([3]) for any choice of G(0). It is in 
fact interesting to observe that if we set 0(0) = a, then the parameter a defines where the 
b{x) that multiplies the jump is evaluated: when a = b{x) is evaluated before the jump, 
while a = 1/2 corresponds to calculating b{x) in the middle of the jump. In the literature 
on GWN, these choices are associated to the I and S prescriptions, respectively [l6|, Il8 |. 
Conversely, as just seen for a discrete jump process, the S interpretation of the SDE ([2]) 
does not correspond to any of the a prescriptions. In other words, there is not an immediate 
intuitive interpretation of the S prescription. 



III. MULTIPLICATIVE COMPOUND POISSON NOISE 



We generalize now our analysis to a process described by the following SDE, 

xit) = a{x,t) + b{x)Cp{^,t), (5) 

where ^p{i^, t) = Yl!i=i wQrit — U) is a colored compound Poisson processes (CP), with jump 
heights w, each time drawn from a generic PDF p{w)^ and {tj} are random times whose 
sequence is drawn from a homogeneous Poisson counting process {N{t),t > 0} of rate u. The 
case in Section 2 corresponds to the special case of a finite deterministic number of jumps. 
As before, the I interpretation consists of taking r = and, should a jump occur at time t, 
of evaluating b{x) at the r.h.s of Eq. ([5]) before the jump occurrence, i.e. x = x{t~), while 
the S interpretation of Eq. ([5]) corresponds to performing the zero limit of the correlation 
time T of the colored Poisson noise. 

The S ME associated with the GLE (jSj) can be derived through the generating function 
of ^^(z^, t) (see Appendix A), or in a more formal way as 

= [ - t) + ,(e-'-'£^'> - PS(,, i), (6) 

where (■) denotes the ensemble average operator. A simpler alternative derivation of the 
ME (Q, can be obtained using the fact that in the S prescription the rules of calculus 
are preserved. Defining the function ri{x) = /^^fry, the ME can be also written as (see 
Appendix B): 

dP''{x,t) dr, , , r pivix) - r]{x')) 



dt 



In the I prescription, x(t) at time t does not depend on the noise ^^~°(z/, t) = ^p(z/, t) at 



the same time 20[ . From this it follows that 

{b{x)CMt)) = {b{x)){U^,t)). (8) 



Therefore, if ([5]) with r = is interpreted in the I sense, we can change the size o: 
from w to b{x)w, and the corresponding ME can be derived without ambiguity 

dP^{x,t) dr. f°° f X - x'\ P\x' ,t) 



the jumps 
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id 



dt 
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Alternatively we achieved a different form of the I ME (Q that is the I analogous of the S 

ME dnD 

where : : is an operator (analogous to the normal order operator in quantum field theory) 
which indicates that all the derivatives must be placed on the left of the expression, i.e. 
: e-'"i^^(^)F(x) := ^^^g ^^(^)"(6(a;)"F(x)). For details see Appendix C. 

When b{x) = b is constant, by using e~''^a^ P^{x,t) = P^{x — bw,t) the I and S MEs 
become coincident, as expected. In Appendix D we also show that, taking the limit — oo, 
(w) — )■ 0, i.e. infinite frequency and infinitesimally small jumps, such that z/(w)^ = D 
remains constant, Eg s. (El) and f lTOj) reduce to the well known I and S Fokker-Planck equation 



(FPE) for GWN [Mm, respectively. 



IV. PRESCRIPTION-INDUCED JUMP DISTRIBUTIONS 

It is clear from the previous MEs ([7]) and Q that the I and S prescriptions of the GLE 

x{t) = a{x,t) + b{x)^p{u,t) (11) 

lead to different MEs. We now want to determine the connection between the two different 
interpretations. Specifically, we seek the two jump PDFs in the I and S interpretation, pi 
and ps, which give rise to the same process. We also seek how to obtain one form when 
the other is given. To this purpose it is sufficient to equate the two MEs, and (E]) for 
simplicity, from which 

As a result, if Eq. ( IT2l) can be solved, given the jumps PDF and choosing the S (I) 
prescription for Eq. (ITTl) . the solutions pi (ps) of Eq. ( |T2l) give the equivalent corresponding 
I (S) GLE and ME. This is one of the main results of the paper and it provides the connection 
between the prescription-induced jump distributions p/ and ps, allowing link the Ito ME 
and the Stratonovich ME corresponding to a GLE with multiplicative white Poisson noise. 

The previous equation however has a solution only when b{x) is a linear function of x. 
To show this we rewrite Eq. (fT2i) as 



Piiy) = J^Psivix) - v{x')) ^ Fix', y), (13) 



where y = {x — x')/b{x'). Because the l.h.s. of Eq. (fT3|) does not depend on x', we must 
have ffr = 0. If the latter condition holds for all ps, then we get b"{x) = whose solution 
is b{x) = kx + c (see Appendix E for details). For other functional shape of b{x) the jumps 
PDF pi{w) depends also on the state of the system, i.e., the dependence on x of pi{w\x) 
cannot be factored out. In this case, is not even clear to what a Stratonovich prescription 
would correspond to. 

Finally we derive the distribution of the impulses that may be measured from the time 
series of the process (see inset in Fig. [2]). In fact, if a random jump (drawn from p{w)) 
occurs at time t, then the size of the impulse that the whole process experiences is yt = 
x{t + dt) — x{t). From the GLE (JSj), we know that with probability udt, x = b{x) w Qr{t)- 
Taking the limit r — )• 0, and using the definition of ri{x) we have (see Appendix B) 

, wb(x), (I) 

yit) = { _ (14) 

7] ^{ri{x) + w) — X (S) 
and thus we obtain 

/+00 -1 
dxdw-—-P'{x, t)p{w)S{w - y/b{x)) (15) 
\Kx)\ 

J-oo \b{x + y)\ 

i.e. the prescriptions characterize the PDF of the impulses of the whole process. 



V. APPLICATION TO SOIL SALINIZATION 

The above mathematical problems naturally arise in the context of the process of soil 
salinization. This is an extremely relevant environmental problem as four million km^ in arid 
and semi-arid lands are affected by soil salinization, causing vegetation dieoff and possible 



desertification 



21 



22| . In natural salinization (unlike the anthropogenic one due to irriga- 
tion), salt may accumulate in surface soils by dry and wet deposition due to wind and rain. 
In this problem, state-dependent Poisson jumps arise naturally when writing the salt mass 
balance equation at the daily-to-monthly time scale for soil root zone used as the control 



volume 



22[. Salt inputs due to rainfall and wind act almost continuously in time, while 



the state-dependent losses of salt occur through negative jumps due to the leaching caused 
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by intense rainfall events. Schematically, the salt mass at time, x{t), in the root zone is 
described by the GLE: 

^ = T-xU^,t), (17) 

where T is the time-averaged salt mass input flux, ^p(z/, t) is the leaching flux toward deeper 
layers, which can be approximated by a WP process with p{w) = fiexp{—fiw)Q{w). The 
leaching parameters u (frequency of leaching events) and n (mean jump) can be expressed 
in terms of the climatic, soil and vegetation properties 22|. Because the typical duration 
of leaching events is on the order of a few hours, while the equilibration times of salt in the 
soil solution (proportional to the inverse of its dissolution rate) tend to be smaller (minutes 
to hours), this means that the inertia in the dynamics is small {a ^ r) and the physically 
correct interpretation is likely to be the Stratonovich one. 
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FIG. 2: (Color online) Comparison of the steady state PDF of Eq. (|17p : S solution (solid line, 
obtained analytically from Eq. (jlSp ). I solution (dash-dot, from numerical simulation), I solution 
using the jump distribution given by Eq. (j2Up (dotted line, from numerical simulation). The 
numerical simulations confirm our analytical results. Inset: simulated trajectory of the salt mass 
under the two different prescriptions. Note that if artificial reflecting barriers are not imposed, 
the salt mass given by the I prescription of Eq. (jl7p may assume unphysical negative values. The 
parameters used for the simulation are \x = 0.463, v = 0.15 day^^ and T = 30mg/day. 
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The stationary solution of Eq. ([6]) in the S prescription is a Gamma distribution (Fig. 

ps^x) = ATe-^'^/^^x^/'^, (18) 

1+;^ 1 0- 

for X > and where Af = ( y) /r(-^) is the normahzation constant and T{z) the complete 
gamma function of argument z. Eq. fllSp summarizes the soil salinity statistics as a function 
of climate, soil and vegetation parameters, which may in turn be used in conjunction with 
the soil moisture statistic to obtain a full characterization of the salt concentration in the 



root zone and the ensuing risk of salinization 221 ]. 



From Eq. f[T6|) it is possible to derive the PDF of the impulses of the process for the S 
interpretation as 

P'iy) = ee^'&i-y), (19) 

which is an exponential distribution controlled by the parameter e = ly/T, given by the ratio 
between the rate of leaching events and the average rate of salt input. Thus if time series 
of the process are available, the Stratonovich assumption can be checked by backtracking 
information on the physical timescales involved in the process, via a comparison with ex- 
perimental data. A further support for the S interpretation of Eq. f|T71) is given by the fact 
that X must remain positive after a jump, a fact that is not ensured by the I interpreta- 
tion unless a reflecting boundary in x = is imposed (see Fig. [2]). We also computed the 
prescription-induced jump distributions correspondence for this case {b{x) = —x), which is 
/95(ln|^|) = |^|p/(l — ^), where x' and x are the variables before and after the jump, 
respectively. By taking into account that in the S prescription x,x' > 0, the I-jump PDF 
equivalent to psiw) = 'ye~^^Q{w) is 

pi{z)=^{l-zy-\ ze[0,l]. (20) 

This equivalence is indeed remarkable because it considerably facilitates the numerical sim- 
ulation of the salinity equation in the S formulation (see Fig. |2]). On the other hand, if 
the GLE ( |T7l) were interpreted in the I sense, the ratio x/x' could also be negative and the 
solution of Eq. ( IT2l) . for pj = 'yQ{w)e~^^ , would read 



Ps{w) = 76 ^ 



eiw)e^''^ + e-^"~" we]-oo, +oo[. (21) 



This implies that possible negative jumps (that occur for x < 0) in the I prescription for the 
given pi{w), would be explicitly present in the corresponding equivalent S-jump PDF ps{w) 
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FIG. 3: (Color online) Comparison between a jump exponential distribution p{w) with mean 
1/7 = 0.8, and the solutions pi{z), ps{z) of the prescription induced jumps corresponding to Eqs. 
(|20l) and (f2T]l respectively, corresponding to the given p{w). 

(see Fig. dS])). 



VI. CONCLUSIONS 

In this work we have proposed a novel approach to solve the Ito-Stratonovich (I-S) 
dilemma for GLE with multiplicative WP noise. We have shown how different interpre- 
tations lead to different results and that choosing between the I and S prescriptions is 
crucial to describe correctly the dynamics of the model systems, and how this choice can be 
determined by physical information about the timescales involved in the process. Moreover, 
we have addressed the related issue of finding a connection between the I and S interpreta- 
tions in the case of linear WP noise. Differently from the introduction of a drift previously 
proposed , we have found such connection in a transformation of the jumps PDFs and 

tested our results numerically. Our results are also consistent with the physics of the random 
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forcing, which takes place at specific points in time, whereas a continuously-acting spurious 
drift would conceptually violate the causality of the process. In particular, once the GLE 
(ITT]) is given, its I and S interpretations are shown to be equivalent if pi and ps satisfy the 
prescription-induced jumps PDF Eq. f[T2]) . The case of nonlinear multiplicative WP noise 
will be studied elsewhere. We have applied our results to the geophysical problem of soil 
salinization, by solving a minimalist model that describes the salt mass and concentration 
in a soil control volume as a function of climatic and ecohydrological parameters. 
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VII. APPENDIX A 



The stochastic process under study is described by the GLE ([5]) presented in the main 
text. For simphcity in the following we have set ^p{t, v) = ^{t). The CP is characterized by 
the correlation structure ( (■) denotes the ensemble average) 



(22) 



where r is the characteristic time of the process and we have omitted all the sub- and super- 
scripts to simplify the notation. If $f is the generating function of CP at time t, then 



(23) 



n=0 

+ 00 

n=0 



-ut 



■ ^ 7 = 1 



Wj Jq v(s)&T{s—tj)ds 



-ut • 



dr I dwp{w) exp 



Wj / v{s)Qt-{s ~ tj)ds 



Moreover if we define p = j _^ ^^^^ p{w)dw as the characteristic function of p(tf), then we 
have 



^t[v\ = exp 



ut + u J drp^J v{s)Qr{s — T)dT^ 



and thus 



"^tlv] = ln^t[v] = u j dr p {^j v{s)Q)r{,s - T)dT ] -1 



(24) 



(25) 



The Stratonovich interpretation of Eq. f lTTl) arises when the limit O^ft — r) — t- 5{t — r) 

nn 

is taken [13|, |17[, that is considering a white Poisson process (WP) as the zero limit of the 
correlation time of the corresponding CP. For a WP the logarithm of the generating function 
thus reads 

\ij^[v] = u dr[p{v{r)) -l]. (26) 
I 

^tH = / dSi---dSnV{Si)---v{s.r,){{^{h)---^{Q))^, (27) 

n=l Jo 



where ((■ ■ ■ ))j is the j-th cumulant, i.e., = (•), ((■■))2 = (••) — etc... 
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From Eqs. (126|) and (1271) we obtain the explicit formula to calculate the cumulants 

mM. = = ^P"(«),„,.^(*. - (29) 

((«<.)■••{(«». = i„fe(i";fo(i„) = ^p""Wl„.4fe-*i)----*.-«»-i).(30) 

In this way, once p{w) is given, we have a complete description of the WP. For example in the 
case of exponential distributed jumps, i.e. p{w) = -j^e "W, the WP is fully characterized 
by the moments 

{m))i = ^{^) (31) 
mtint2)))2 = v{w^)5{u-h) (32) 

mh)-^{tn)))n = iy{w^)6{U - h)...6{tn-l - tn)- (33) 

Once we have calculate all the moments of the WP process we can easily achieve the ME 
corresponding to the GLE (fTTl) . For a given realization of ^ the solution of Eq. ([5]) is 

/(x,t|0 = 5(a;-x(t)). (34) 

To obtain the general solution of Eq. ffTTl) we simply have to take the ensemble average of 
different trajectories 

(/(x,t|0)=P''(a;,t). (35) 
Differentiating both sides of Eq. fl34l) and using Eq. ffTTj) we have 

aip^(x,t|0 = dj{x - x{t))[-x{t)] (36) 
= -dj{x - x{t))[a{x{t),t) + b{x{t))i{t)] (37) 
= -dj{x - x{t))[a{x{t),t) + h{x{t))i{t)], (38) 

and thus we obtain a forward ME for the PDF conditioned by a given realization of the WP 

= -0{x,d,,t)p'{x,m, (39) 

where 0{x,dx,t) = dx[a{x(t),t) + b{x{t))$,(t)] is the forward time evolution operator. The 

solution of Eq. f l5I?]) . for the initial condition p^{x{0), 0|^) = 6{x — a;(0)) is 

ft 



p'{x,t\0 = T(exp[ - / {d,a{x,r) + dXx{rmr))dT] ]6{x-x{0)), (40) 

14 



where T is the T-product operator. Using Eq. (135|) and the Kubo relation (127|) . an exphcit 
formula for the general formal solution of the GLE ffTTj) in the Stratonovich prescription is 
obtained 



exp [ — / dxa{x, T)dT — / dti - ■ ■ dt^dx 
Jo Jo Jo 



6(x(ti))---9,6(x(t„)) X (41) 

(42) 

Thanks to Eqs. 0281) . 0291) and 03Up we have a complete characterization of the cumulants, 
and thus substituting Eq. (I5U]) into Eq. (HTj) we obtain 

P'^ix.t) = T^exp[-^ dxa{x,r)dT + fy (-5,.6(x(r)))" p(0)^"^rfr] 

= T (^exp [ - l^dMx,r)dT - £ rfr(e-^^''(^«) - 1>^(^J^ (43) 



Eventually, differentiating Eq. f l43|) with respect to t we obtain the ME corresponding to 
the GLE (fTT!) in the Stratonovich interpretation 



di 



d_ 

dx 



j(x,t) + z/(e-'"^'(") -1) 



p{w) 



P^(x,t), 



(44) 



that is the ME reported in the main text. 



VIII. APPENDIX B 

We now show the derivation of the S ME ([7]) in the main text and its equivalence with 
Eq. (E]). We can write the GLE ([5]) as 

a(x,t), with probability 1 — u dt] 

x{t) = { (45) 

b{x) w hT-{t), with probability u dt, 



where hrit) = Qrit)- We now consider only the effect of the jumps on x. From Eq. (H5|) 
we have that dx/b{x{t)) = whT-{t)dt, and setting 

^'^'^ - ' - Vix) = r (46) 



dx b{x) J b{x') 

Eq. f H5l) becomes 

drj{x{t)) = w h^-dt, (47) 
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which integrated between t and t + dt reads 

r]{x{t + dt)) = r/(x(t)) + w AQ^it) x{t + dt) = 7]'^[r]{x{t)) + wAe^(t)], (48) 

where AQ^{t) = Q^{t + dt) - Q^{t). 

Finally we can write the discrete ME corresponding to the GLE ffTTl) interpreted in the 
Stratonovich sense 

POO 

P^{x,t + dt) = {l-v)dt I dx' P^{x\t)5{x -{a{x')dt + x)) + (49) 

Jo 



oo roo 



+ udt / / p{w)P^{x',t)6 {x - [t] ^[vix') + w])) dwdx', 
Jo Jo 

where we have performed the limit r — )■ of the GLE ([5]) and used the fact that 
\imr-^o AQrit) = 1. The integral in the r.h.s of Eq. ( l49l) can be rewritten, invert- 
ing the Dirac Delta with respect to w and using the rule of the inverse function, as 
Io° lo p{'w)P^{x') dw dx' and thus, after taking the continuum time limit, 

the Master Equation fj49|) becomes 

= Hx,t)P'{x^t)] + uj^ ^^^^^^~^J^^^^^ P^(x-,t)c^x- - .P^(x,t), (50) 

that is Eq. (JTj) reported in the main text. 

In order to show the equivalence between Eqs. ([6]) and (JTj) we define 

Q{x,w) = / P''^{x')6 {x- {7]-^[r]{x') + w])) dx', (51) 



so we have that the integral in Eq. ( !50|) is simply Q{x,w)p{w)dw. 

Differentiating Eq. ( ISTl) with respect to w we obtain the partial differential equation 

dwQ{x,w) = -dxb{x)Q{x,w) = -T-LQ, (52) 

where we used Eqs. ( H6|) and (IST!) and the definition of the derivative of the inverse function. 
The solution of Eq. ( !52l) is 

Q{x,w) = e-'"^Q{x,0) = e-"'^-''(^)p^(x). (53) 

We thus have 

PM^l^I^^^pS^^'^ t)dx' = {e~-9-'^^^)P'{x, t), (54) 

which substituted in Eq. (!50|) proves the equivalence between the MEs ([6]) and ([7j) in the 
main text. 
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IX. APPENDIX C 



We present in this Appendix the derivation for the I ME (fTOj) and its equivalence with 
Eq. We first note that the integral in the ME can be rewritten as 



p\ — — — — -— — '—dx'= I I p(w)6(x — x' — wb(x'))P (x' ,t)dx'dw (55) 

, V K^ ) ) J-ooJ-oo 

Formally expanding the Dirac delta 

6(x -x' - wb(x')) = V ^^^(-^rb(xTHx - x') (56) 

n\ ox 

n=0 

and substituting Eq. (ISBj) in (13^ we have 

where {■.e-'"i^^^^^:)P\x,t) = E^^o ^^(^)"^(a;)"i^^(a;, t). Using the expression §7i) in 
Eq. (El) we obtain the ME 

X. APPENDIX D 

We now derive the well known FPE corresponding to the GLE ([5]) when ^(t) is a GWN 
with mean {^{t)) = and correlation {^{t) ^{s)) = 2D 5{t - s), from the MEs ^ and ^ 
presented in the main text. We generalize our results to any jump size PDF of the form 

p[w) = 7/(7^^), (58) 

with 7 > and J w'"'p{w)dw = {w^)p < 00 Vn. We note that the latter condition implies 

7 / dww'^f{'yw) = 7""/ dzz'^f{z) = 7""(2")/ < 00 Vn. 

Stratonovich Eq. The case for the Stratonovich prescription has been first presented in 
[3] . The FPE corresponding to multiplicative GWN process interpreted in the Stratonovich 
sense is 

^P^(x,t) = -^[a{x,t)P'{x,t)] + D-^b{x)^b{x)P'{x,t). (59) 

Once we consider a zero mean WP process, the ME ([6]) reads as [tI 
dP^{x,t) d 



dt dx 



t) - u{w)b{x)\P^{x,t) 



+ z/(e-"'a^''W - 1) p^(x,t) (60) 



.^K^,t)p^(x,t)] + (61) 

n=l 
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where the integral in the r.h.s of Eq. (160|1 has been expanded as 

n=0 n=0 ' 

(62) 

Taking the hmit i/, 7 — i- oo, such that ^ = D', then ^ foi n > 2 and the latter ME 
fiGTll corresponds exactly to the FPE (JSl) with D = D'^^. 

ltd Eq. The FPE corresponding to multiplicative GWN process interpreted with the Ito 
prescription is 

^P^(x, t) = [a{x, t) P{x, t)]+D-^ [b{xY P\x, t)] . (63) 

We now repeat the same procedure as before, starting from the zero mean I ME 

^^^^ = [H^,t) - u{w)bix)]P\x,t)] + e-^''^-) :) - l)P\x,t). (64) 

We can expand the r.h.s. remembering that the operator :: means that all the derivatives 
must be placed on the left of the expression 

e-^^(^) ■.)-l)P\x,t) = -^^(i/^)p^[&(^)P'(x,t)]+z.X](--)"^(5^)"[^(^)"^'(^'^)] 



n=2 



(65) 



Eventually, inserting Eq. ( 165|) in the I ME ( l64l) and taking i/, 7 — 00 with = D' and 
D = D'^^ we obtain the I FPE (jeS]). 



Appendix E 



In this section we show how a solution of Eq. f|T2|) . rewritten as 



Pi{y) = ^-J^Psivix) - v{x')) = F{x\ y), (66) 



where y = [x — x')/b{x'), exist only if is 6 is a linear function. 

dx' 



Because the l.h.s. of Eq. (!66|) does not depend on x', we must have |^ = 0, that explicitly 



read as 



= psivix)-v{x'))[sgn[b{x')]^^^-sgn[b{x)]^^ 

+ p'5(r](x)-r/(x'))i^[r/'(x)(l + y6'(x)-r^V))]- (67) 
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The latter, using Eq. (HGj) . can be expressed as 

Eq. (!68l) must hold for all p5, then the solution of Eq. ( |68|) is given by the function 6 
that satisfies the conditions 

b{x'){l+yb'{x')) = b{x) (69) 
b{x')b'{x){l + yb'ix')) = b{x)b'{x'). (70) 

Combining Eqs. fl69|) and flTOj) and using a; = b{x')y + x' we obtain the equation h'{b{x')y + 
x') = b'{x'). If we take the derivative of both side with respect to the independent variable y, 
then we have b" {b{x')y + x')b{x') = 0. This implies b"{x) = Vx, which solution is b{x) = kx 
(with k any constant). 
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